Graph Perceiver Network for Lung Tumor and Bronchial Premalignant Lesion Stratification from Histopathology.

Bronchial premalignant lesions (PMLs) precede the development of invasive lung squamous cell carcinoma (LUSC), posing a significant challenge in distinguishing those likely to advance to LUSC from those that might regress without intervention. In this context, we present a novel computational approach, the Graph Perceiver Network, leveraging hematoxylin and eosin-stained whole slide images to stratify endobronchial biopsies of PMLs across a spectrum from normal to tumor lung tissues. The Graph Perceiver Network outperforms existing frameworks in classification accuracy predicting LUSC, lung adenocarcinoma, and nontumor (normal) lung tissue on The Cancer Genome Atlas and Clinical Proteomic Tumor Analysis Consortium datasets containing lung resection tissues while efficiently generating pathologist-aligned, class-specific heat maps. The network was further tested using endobronchial biopsies from two data cohorts, containing normal to carcinoma in situ histology, and it demonstrated a unique capability to differentiate carcinoma in situ lung squamous PMLs based on their progression status to invasive carcinoma. The network may have utility in stratifying PMLs for chemoprevention trials or more aggressive follow-up.

Graph attention-based fusion of pathology images and gene expression for prediction of cancer survival.

Multimodal machine learning models are being developed to analyze pathology images and other modalities, such as gene expression, to gain clinical and biological insights. However, most frameworks for multimodal data fusion do not fully account for the interactions between different modalities. Here, we present an attention-based fusion architecture that integrates a graph representation of pathology images with gene expression data and concomitantly learns from the fused information to predict patient-specific survival. In our approach, pathology images are represented as undirected graphs, and their embeddings are combined with embeddings of gene expression signatures using an attention mechanism to stratify tumors by patient survival. We show that our framework improves the survival prediction of human non-small cell lung cancers, outperforming existing state-of-the-art approaches that leverage multimodal data. Our framework can facilitate spatial molecular profiling to identify tumor heterogeneity using pathology images and gene expression data, complementing results obtained from more expensive spatial transcriptomic and proteomic technologies.

Integrative genetic and genomic networks identify microRNA associated with COPD and ILD.

Chronic obstructive pulmonary disease (COPD) and interstitial lung disease (ILD) are clinically and molecularly heterogeneous diseases. We utilized clustering and integrative network analyses to elucidate roles for microRNAs (miRNAs) and miRNA isoforms (isomiRs) in COPD and ILD pathogenesis. Short RNA sequencing was performed on 351 lung tissue samples of COPD (n = 145), ILD (n = 144) and controls (n = 64). Five distinct subclusters of samples were identified including 1 COPD-predominant cluster and 2 ILD-predominant clusters which associated with different clinical measurements of disease severity. Utilizing 262 samples with gene expression and SNP microarrays, we built disease-specific genetic and expression networks to predict key miRNA regulators of gene expression. Members of miR-449/34 family, known to promote airway differentiation by repressing the Notch pathway, were among the top connected miRNAs in both COPD and ILD networks. Genes associated with miR-449/34 members in the disease networks were enriched among genes that increase in expression with airway differentiation at an air-liquid interface. A highly expressed isomiR containing a novel seed sequence was identified at the miR-34c-5p locus. 47% of the anticorrelated predicted targets for this isomiR were distinct from the canonical seed sequence for miR-34c-5p. Overexpression of the canonical miR-34c-5p and the miR-34c-5p isomiR with an alternative seed sequence down-regulated NOTCH1 and NOTCH4. However, only overexpression of the isomiR down-regulated genes involved in Ras signaling such as CRKL and GRB2. Overall, these findings elucidate molecular heterogeneity inherent across COPD and ILD patients and further suggest roles for miR-34c in regulating disease-associated gene-expression.

Single-Cell Characterization of Pulmonary Nodules Implicates Suppression of Immunosurveillance across Early Stages of Lung Adenocarcinoma.

A greater understanding of molecular, cellular, and immunological changes during the early stages of lung adenocarcinoma development could improve diagnostic and therapeutic approaches in patients with pulmonary nodules at risk for lung cancer. To elucidate the immunopathogenesis of early lung tumorigenesis, we evaluated surgically resected pulmonary nodules representing the spectrum of early lung adenocarcinoma as well as associated normal lung tissues using single-cell RNA sequencing and validated the results by flow cytometry and multiplex immunofluorescence (MIF). Single-cell transcriptomics revealed a significant decrease in gene expression associated with cytolytic activities of tumor-infiltrating natural killer and natural killer T cells. This was accompanied by a reduction in effector T cells and an increase of CD4+ regulatory T cells (Treg) in subsolid nodules. An independent set of resected pulmonary nodules consisting of both adenocarcinomas and associated premalignant lesions corroborated the early increment of Tregs in premalignant lesions compared with the associated normal lung tissues by MIF. Gene expression analysis indicated that cancer-associated alveolar type 2 cells and fibroblasts may contribute to the deregulation of the extracellular matrix, potentially affecting immune infiltration in subsolid nodules through ligand-receptor interactions. These findings suggest that there is a suppression of immune surveillance across the spectrum of early-stage lung adenocarcinoma. SIGNIFICANCE: Analysis of a spectrum of subsolid pulmonary nodules by single-cell RNA sequencing provides insights into the immune regulation and cell-cell interactions in the tumor microenvironment during early lung tumor development.

Convergence of YAP/TAZ, TEAD and TP63 activity is associated with bronchial premalignant severity and progression.

BACKGROUND: Bronchial premalignant lesions (PMLs) are composed primarily of cells resembling basal epithelial cells of the airways, which through poorly understood mechanisms have the potential to progress to lung squamous cell carcinoma (LUSC). Despite ongoing efforts that have mapped gene expression and cell diversity across bronchial PML pathologies, signaling and transcriptional events driving malignancy are poorly understood. Evidence has suggested key roles for the Hippo pathway effectors YAP and TAZ and associated TEAD and TP63 transcription factor families in bronchial basal cell biology and LUSC. In this study we examine the functional association of YAP/TAZ, TEADs and TP63 in bronchial epithelial cells and PMLs. METHODS: We performed RNA-seq in primary human bronchial epithelial cells following small interfering RNA (siRNA)-mediated depletion of YAP/TAZ, TEADs or TP63, and combined these data with ChIP-seq analysis of these factors. Directly activated or repressed genes were identified and overlapping genes were profiled across gene expression data obtained from progressive or regressive human PMLs and across lung single cell RNA-seq data sets. RESULTS: Analysis of genes regulated by YAP/TAZ, TEADs, and TP63 in human bronchial epithelial cells revealed a converged transcriptional network that is strongly associated with the pathological progression of bronchial PMLs. Our observations suggest that YAP/TAZ-TEAD-TP63 associate to cooperatively promote basal epithelial cell proliferation and repress signals associated with interferon responses and immune cell communication. Directly repressed targets we identified include the MHC Class II transactivator CIITA, which is repressed in progressive PMLs and associates with adaptive immune responses in the lung. Our findings provide molecular insight into the control of gene expression events driving PML progression, including those contributing to immune evasion, offering potential new avenues for lung cancer interception. CONCLUSIONS: Our study identifies important gene regulatory functions for YAP/TAZ-TEAD-TP63 in the early stages of lung cancer development, which notably includes immune-suppressive roles, and suggest that an assessment of the activity of this transcriptional complex may offer a means to identify immune evasive bronchial PMLs and serve as a potential therapeutic target.

Interrogating the precancerous evolution of pathway dysfunction in lung squamous cell carcinoma using XTABLE.

Lung squamous cell carcinoma (LUSC) is a type of lung cancer with a dismal prognosis that lacks adequate therapies and actionable targets. This disease is characterized by a sequence of low- and high-grade preinvasive stages with increasing probability of malignant progression. Increasing our knowledge about the biology of these premalignant lesions (PMLs) is necessary to design new methods of early detection and prevention, and to identify the molecular processes that are key for malignant progression. To facilitate this research, we have designed XTABLE (Exploring Transcriptomes of Bronchial Lesions), an open-source application that integrates the most extensive transcriptomic databases of PMLs published so far. With this tool, users can stratify samples using multiple parameters and interrogate PML biology in multiple manners, such as two- and multiple-group comparisons, interrogation of genes of interests, and transcriptional signatures. Using XTABLE, we have carried out a comparative study of the potential role of chromosomal instability scores as biomarkers of PML progression and mapped the onset of the most relevant LUSC pathways to the sequence of LUSC developmental stages. XTABLE will critically facilitate new research for the identification of early detection biomarkers and acquire a better understanding of the LUSC precancerous stages.

Lung squamous cell carcinoma is the second most common lung cancer. However, very little is known about how normal tissues in the lung develop in to these tumours. Like many cancers, this transformation comprises of an intermediate phase where healthy cells begin to form lesions that may (or may not) progress in to tumours. Understanding the biology of these lesions in lung squamous cell carcinoma may help clinicians detect them before they become cancerous. Knowing which genes are switched on and off during this intermediary phase can provide clues as to how these lesions form. There are already some publicly available transcriptional datasets showing the activity of tens of thousands of genes in pre-cancerous lesions extracted from patients with lung squamous cell carcinoma. But not every laboratory has the bioinformatic tools and skills required to interrogate these extensive databases. To address this, Roberts et al. built an open-source platform called XTABLE (short for Exploring Transcriptomes of Bronchial Lesions) which can analyse transcriptional datasets in multiple ways depending on the needs of the user. For instance, the tool can stratify the data into groups based on different parameters, such as the lesions potential to progress in to cancer, to see how the genes of the groups compare. It can also analyse the activity of individual genes and sets of genes involved in the same biological processes. Using XTABLE, Roberts et al. showed that a biological process linked to lung squamous cell carcinoma is also involved in the formation of pre-cancerous lesions. This suggests that molecules and genes associated with this process could potentially help scientists design prevention strategies. XTABLE will help researchers to better understand the biology of pre-cancerous lesions and how they develop in to tumours. Moreover, it will make it easier for scientists to validate their hypotheses using data collected from patients. The tool could also be useful for scientists interested in other types of lung cancers that share a similar biology.

Inactivation of LATS1/2 drives luminal-basal plasticity to initiate basal-like mammary carcinomas.

Basal-like breast cancers, an aggressive breast cancer subtype that has poor treatment options, are thought to arise from luminal mammary epithelial cells that undergo basal plasticity through poorly understood mechanisms. Using genetic mouse models and ex vivo primary organoid cultures, we show that conditional co-deletion of the LATS1 and LATS2 kinases, key effectors of Hippo pathway signaling, in mature mammary luminal epithelial cells promotes the development of Krt14 and Sox9-expressing basal-like carcinomas that metastasize over time. Genetic co-deletion experiments revealed that phenotypes resulting from the loss of LATS1/2 activity are dependent on the transcriptional regulators YAP/TAZ. Gene expression analyses of LATS1/2-deleted mammary epithelial cells notably revealed a transcriptional program that associates with human basal-like breast cancers. Our study demonstrates in vivo roles for the LATS1/2 kinases in mammary epithelial homeostasis and luminal-basal fate control and implicates signaling networks induced upon the loss of LATS1/2 activity in the development of basal-like breast cancer.

Smoking modulates different secretory subpopulations expressing SARS-CoV-2 entry genes in the nasal and bronchial airways.

SARS-CoV-2 infection and disease severity are influenced by viral entry (VE) gene expression patterns in the airway epithelium. The similarities and differences of VE gene expression (ACE2, TMPRSS2, and CTSL) across nasal and bronchial compartments have not been fully characterized using matched samples from large cohorts. Gene expression data from 793 nasal and 1673 bronchial brushes obtained from individuals participating in lung cancer screening or diagnostic workup revealed that smoking status (current versus former) was the only clinical factor significantly and reproducibly associated with VE gene expression. The expression of ACE2 and TMPRSS2 was higher in smokers in the bronchus but not in the nose. scRNA-seq of nasal brushings indicated that ACE2 co-expressed genes were highly expressed in club and C15orf48+ secretory cells while TMPRSS2 co-expressed genes were highly expressed in keratinizing epithelial cells. In contrast, these ACE2 and TMPRSS2 modules were highly expressed in goblet cells in scRNA-seq from bronchial brushings. Cell-type deconvolution of the gene expression data confirmed that smoking increased the abundance of several secretory cell populations in the bronchus, but only goblet cells in the nose. The association of ACE2 and TMPRSS2 with smoking in the bronchus is due to their high expression in goblet cells which increase in abundance in current smoker airways. In contrast, in the nose, these genes are not predominantly expressed in cell populations modulated by smoking. In individuals with elevated lung cancer risk, smoking-induced VE gene expression changes in the nose likely have minimal impact on SARS-CoV-2 infection, but in the bronchus, smoking may lead to higher viral loads and more severe disease.

Clinical Study of Aspirin and Zileuton on Biomarkers of Tobacco-Related Carcinogenesis in Current Smokers.

The chemopreventive effect of aspirin and other non-steroidal anti-inflammatory drugs (NSAIDs) on lung cancer risk is supported by epidemiologic and preclinical studies. Zileuton, a 5-lipoxygenase inhibitor, has additive activity with NSAIDs against tobacco carcinogenesis in preclinical models. We hypothesized that cyclooxygenase plus 5-lipoxygenase inhibition would be more effective than a placebo in modulating the nasal epithelium gene signatures of tobacco exposure and lung cancer. We conducted a randomized, double-blinded study of low-dose aspirin plus zileuton vs. double placebo in current smokers to compare the modulating effects on nasal gene expression and arachidonic acid metabolism. In total, 63 participants took aspirin 81 mg daily plus zileuton (Zyflo CR) 600 mg BID or the placebo for 12 weeks. Nasal brushes from the baseline, end-of-intervention, and one-week post intervention were profiled via microarray. Aspirin plus zilueton had minimal effects on the modulation of the nasal or bronchial gene expression signatures of smoking, lung cancer, and COPD but favorably modulated a bronchial gene expression signature of squamous dysplasia. Aspirin plus zileuton suppressed urinary leukotriene but not prostaglandin E2, suggesting shunting through the cyclooxygenase pathway when combined with 5-lipoxygenase inhibition. Continued investigation of leukotriene inhibitors is needed to confirm these findings, understand the long-term effects on the airway epithelium, and identify the safest, optimally dosed agents.

A Graph-Transformer for Whole Slide Image Classification.

Deep learning is a powerful tool for whole slide image (WSI) analysis. Typically, when performing supervised deep learning, a WSI is divided into small patches, trained and the outcomes are aggregated to estimate disease grade. However, patch-based methods introduce label noise during training by assuming that each patch is independent with the same label as the WSI and neglect overall WSI-level information that is significant in disease grading. Here we present a Graph-Transformer (GT) that fuses a graph-based representation of an WSI and a vision transformer for processing pathology images, called GTP, to predict disease grade. We selected 4,818 WSIs from the Clinical Proteomic Tumor Analysis Consortium (CPTAC), the National Lung Screening Trial (NLST), and The Cancer Genome Atlas (TCGA), and used GTP to distinguish adenocarcinoma (LUAD) and squamous cell carcinoma (LSCC) from adjacent non-cancerous tissue (normal). First, using NLST data, we developed a contrastive learning framework to generate a feature extractor. This allowed us to compute feature vectors of individual WSI patches, which were used to represent the nodes of the graph followed by construction of the GTP framework. Our model trained on the CPTAC data achieved consistently high performance on three-label classification (normal versus LUAD versus LSCC: mean accuracy = 91.2 ± 2.5%) based on five-fold cross-validation, and mean accuracy = 82.3 ± 1.0% on external test data (TCGA). We also introduced a graph-based saliency mapping technique, called GraphCAM, that can identify regions that are highly associated with the class label. Our findings demonstrate GTP as an interpretable and effective deep learning framework for WSI-level classification.

Pathologic and gene expression comparison of CT- screen detected and routinely detected stage I/0 lung adenocarcinoma in NCCN risk-matched cohorts.

INTRODUCTION: Although three randomized control trials have proven mortality benefit of CT lung cancer screening (CTLS), <5% of eligible US smokers are screened. Some attribute this to fear of harm conveyed at shared decision visits, including the harm of overdiagnosis/overtreatment of indolent BAC-like adenocarcinoma. METHODS: Since the frequency of indolent cancers has not been compared between CTLS and routinely detected cohorts, we compare pathology and RNA expression of 86 NCCN high-risk CTLS subjects to 83 high-risk (HR-R) and 51 low-risk (LR-R) routinely detected patients. Indolent adenocarcinoma was defined as previously described for low malignant potential (LMP) adenocarcinoma along with AIS/MIA. Exome RNA sequencing was performed on a subset of high-risk (CTLS and HR-R) FFPE tumor samples. RESULTS: Indolent adenocarcinoma (AIS, MIA, and LMP) showed 100% disease-specific survival (DSS) with similar frequency in CTLS (18%) and HR-R (20%) which were comparatively lower than LR-R (33%). Despite this observation, CTLS exhibited intermediate DSS between HR-R and LR-R (5-year DSS: 88% CTLS, 82% HR-R, & 95% LR-R, p = 0.047), possibly reflecting a 0.4 cm smaller median tumor size and lower frequency of tumor necrosis compared to HR-R. WGCNA gene modules derived from TCGA lung adenocarcinoma correlated with aggressive histologic patterns, mitotic activity, and tumor invasive features, but no significant differential expression between CTLS and HR-R was observed. CONCLUSION: CTLS subjects are at no greater risk of overdiagnosis from indolent adenocarcinoma (AIS, MIA, and LMP) than risk-matched patients whose cancers are discovered in routine clinical practice. Improved outcomes likely reflect detection and treatment at smaller size.

Smoking Modulates Different Secretory Subpopulations Expressing SARS-CoV-2 Entry Genes in the Nasal and Bronchial Airways.

Background : SARS-CoV-2 infection and disease severity are influenced by viral entry (VE) gene expression patterns in airway epithelium. The similarities and differences of VE gene expression (ACE2, TMPRSS2, and CTSL) across nasal and bronchial compartments has not been fully characterized using matched samples from large cohorts. Results : Gene expression data from 793 nasal and 1,673 bronchial brushes obtained from individuals participating in lung cancer screening or diagnostic workup revealed that smoking was the only clinical factor significantly and reproducibly associated with VE gene expression. ACE2 and TMPRSS2 expression were higher in smokers in the bronchus but not in the nose. scRNA-seq of nasal brushings indicated that ACE2 co-expressed genes were highly expressed in club and C15orf48 + secretory cells while TMPRSS2 co-expressed genes were highly expressed in keratinizing epithelial cells. In contrast, these ACE2 and TMPRSS2 modules were highly expressed in goblet cells in scRNA-seq from bronchial brushings. Cell-type deconvolution of the RNA-seq confirmed that smoking increased the abundance of several secretory cell populations in the bronchus, but only goblet cells in the nose. Conclusions : The association of ACE2 and TMPRSS2 with smoking in the bronchus is due to their high expression in goblet cells which increase in abundance in current smoker airways. In contrast, in the nose these genes are not predominantly expressed in cell populations modulated by smoking. Smoking-induced VE gene expression changes in the nose likely has minimal impact on SARS-CoV-2 infection, but in the bronchus, smoking may lead to higher viral loads and more severe disease.

Elevated T cell repertoire diversity is associated with progression of lung squamous cell premalignant lesions.

OBJECTIVE: The immune response to invasive carcinoma has been the focus of published work, but little is known about the adaptive immune response to bronchial premalignant lesions (PMLs), precursors of lung squamous cell carcinoma. This study was designed to characterize the T cell receptor (TCR) repertoire in PMLs and its association with clinical, pathological, and molecular features. METHODS: Endobronchial biopsies (n=295) and brushings (n=137) from high-risk subjects (n=50), undergoing lung cancer screening at approximately 1-year intervals via autofluorescence bronchoscopy and CT, were profiled by RNA-seq. We applied the TCR Repertoire Utilities for Solid Tissue/Tumor tool to the RNA-seq data to identify TCR CDR3 sequences across all samples. In the biopsies, we measured the correlation of TCR diversity with previously derived immune-associated PML transcriptional signatures and PML outcome. We also quantified the spatial and temporal distribution of shared and clonally expanded TCRs. Using the biopsies and brushes, the ratio of private (ie, found in one patient only) and public (ie, found in two or more patients) TCRs was quantified, and the CDR3 sequences were compared with those found in curated databases with known antigen specificities. RESULTS: We detected 39,303 unique TCR sequences across all samples. In PML biopsies, TCR diversity was negatively associated with a transcriptional signature of T cell mediated immune activation (p=4e-4) associated with PML outcome. Additionally, in lesions of the proliferative molecular subtype, TCR diversity was decreased in regressive versus progressive/persistent PMLs (p=0.045). Within each patient, TCRs were more likely to be shared between biopsies sampled at the same timepoint than biopsies sampled at the same anatomic location at different times. Clonally expanded TCRs, within a biopsied lesion, were more likely to be expanded at future time points than non-expanded clones. The majority of TCR sequences were found in a single sample, with only 3396 (8.6%) found in more than one sample and 1057 (2.7%) found in two or more patients (ie, public); however, when compared with a public database of CDR3 sequences, 4543 (11.6%) of TCRs were identified as public. TCRs with known antigen specificities were enriched among public TCRs (p<0.001). CONCLUSIONS: Decreased TCR diversity may reflect nascent immune responses that contribute to PML elimination. Further studies are needed to explore the potential for immunoprevention of PMLs.

Yap/Taz inhibit goblet cell fate to maintain lung epithelial homeostasis.

Proper lung function relies on the precise balance of specialized epithelial cells that coordinate to maintain homeostasis. Herein, we describe essential roles for the transcriptional regulators YAP/TAZ in maintaining lung epithelial homeostasis, reporting that conditional deletion of Yap and Wwtr1/Taz in the lung epithelium of adult mice results in severe defects, including alveolar disorganization and the development of airway mucin hypersecretion. Through in vivo lineage tracing and in vitro molecular experiments, we reveal that reduced YAP/TAZ activity promotes intrinsic goblet transdifferentiation of secretory airway epithelial cells. Global gene expression and chromatin immunoprecipitation sequencing (ChIP-seq) analyses suggest that YAP/TAZ act cooperatively with TEA domain (TEAD) transcription factors and the NuRD complex to suppress the goblet cell fate program, directly repressing the SPDEF gene. Collectively, our study identifies YAP/TAZ as critical factors in lung epithelial homeostasis and offers molecular insight into the mechanisms promoting goblet cell differentiation, which is a hallmark of many lung diseases.

Aberrant epithelial polarity cues drive the development of precancerous airway lesions.

Molecular events that drive the development of precancerous lesions in the bronchial epithelium, which are precursors of lung squamous cell carcinoma (LUSC), are poorly understood. We demonstrate that disruption of epithelial cellular polarity, via the conditional deletion of the apical determinant Crumbs3 (Crb3), initiates and sustains precancerous airway pathology. The loss of Crb3 in adult luminal airway epithelium promotes the uncontrolled activation of the transcriptional regulators YAP and TAZ, which stimulate intrinsic signals that promote epithelial cell plasticity and paracrine signals that induce basal-like cell growth. We show that aberrant polarity and YAP/TAZ-regulated gene expression associates with human bronchial precancer pathology and disease progression. Analyses of YAP/TAZ-regulated genes further identified the ERBB receptor ligand Neuregulin-1 (NRG1) as a key transcriptional target and therapeutic targeting of ERBB receptors as a means of preventing and treating precancerous cell growth. Our observations offer important molecular insight into the etiology of LUSC and provides directions for potential interception strategies of lung cancer.

Characterizing smoking-induced transcriptional heterogeneity in the human bronchial epithelium at single-cell resolution.

The human bronchial epithelium is composed of multiple distinct cell types that cooperate to defend against environmental insults. While studies have shown that smoking alters bronchial epithelial function and morphology, its precise effects on specific cell types and overall tissue composition are unclear. We used single-cell RNA sequencing to profile bronchial epithelial cells from six never and six current smokers. Unsupervised analyses led to the characterization of a set of toxin metabolism genes that localized to smoker ciliated cells, tissue remodeling associated with a loss of club cells and extensive goblet cell hyperplasia, and a previously unidentified peri-goblet epithelial subpopulation in smokers who expressed a marker of bronchial premalignant lesions. Our data demonstrate that smoke exposure drives a complex landscape of cellular alterations that may prime the human bronchial epithelium for disease.

Effect of Intermittent Versus Continuous Low-Dose Aspirin on Nasal Epithelium Gene Expression in Current Smokers: A Randomized, Double-Blinded Trial.

A chemopreventive effect of aspirin (ASA) on lung cancer risk is supported by epidemiologic and preclinical studies. We conducted a randomized, double-blinded study in current heavy smokers to compare modulating effects of intermittent versus continuous low-dose ASA on nasal epithelium gene expression and arachidonic acid (ARA) metabolism. Fifty-four participants were randomized to intermittent (ASA 81 mg daily for one week/placebo for one week) or continuous (ASA 81 mg daily) for 12 weeks. Low-dose ASA suppressed urinary prostaglandin E2 metabolite (PGEM; change of -4.55 ± 11.52 from baseline 15.44 ± 13.79 ng/mg creatinine for arms combined, P = 0.02), a surrogate of COX-mediated ARA metabolism, but had minimal effects on nasal gene expression of nasal or bronchial gene-expression signatures associated with smoking, lung cancer, and chronic obstructive pulmonary disease. Suppression of urinary PGEM correlated with favorable changes in a smoking-associated gene signature (P < 0.01). Gene set enrichment analysis (GSEA) showed that ASA intervention led to 1,079 enriched gene sets from the Canonical Pathways within the Molecular Signatures Database. In conclusion, low-dose ASA had minimal effects on known carcinogenesis gene signatures in nasal epithelium of current smokers but results in wide-ranging genomic changes in the nasal epithelium, demonstrating utility of nasal brushings as a surrogate to measure gene-expression responses to chemoprevention. PGEM may serve as a marker for smoking-associated gene-expression changes and systemic inflammation. Future studies should focus on NSAIDs or agent combinations with broader inhibition of pro-inflammatory ARA metabolism to shift gene signatures in an anti-carcinogenic direction.

Immune evasion before tumour invasion in early lung squamous carcinogenesis.

Early detection and treatment are critical for improving the outcome of patients with cancer1. Understanding the largely uncharted biology of carcinogenesis requires deciphering molecular processes in premalignant lesions, and revealing the determinants of the intralesional immune reaction during cancer development. The adaptive immune response within tumours has previously been shown to be strongest at the earliest stage of carcinoma2,3. Here we show that immune activation and immune escape occur before tumour invasion, and reveal the relevant immune biomarkers of the pre-invasive stages of carcinogenesis in the lung. We used gene-expression profiling and multispectral imaging to analyse a dataset of 9 morphological stages of the development of lung squamous cell carcinoma, which includes 122 well-annotated biopsies from 77 patients. We identified evolutionary trajectories of cancer and immune pathways that comprise (1) a linear increase in proliferation and DNA repair from normal to cancerous tissue; (2) a transitory increase of metabolism and early immune sensing, through the activation of resident immune cells, in low-grade pre-invasive lesions; (3) the activation of immune responses and immune escape through immune checkpoints and suppressive interleukins from high-grade pre-invasive lesions; and, ultimately, (4) the activation of the epithelial-mesenchymal transition in the invasive stage of cancer. We propose that carcinogenesis in the lung involves a dynamic co-evolution of pre-invasive bronchial cells and the immune response. These findings highlight the need to develop immune biomarkers for early detection as well as immunotherapy-based chemopreventive approaches for individuals who are at high risk of developing lung cancer.

Tobacco-Related Alterations in Airway Gene Expression are Rapidly Reversed Within Weeks Following Smoking-Cessation.

The physiologic response to tobacco smoke can be measured by gene-expression profiling of the airway epithelium. Temporal resolution of kinetics of gene-expression alterations upon smoking-cessation might delineate distinct biological processes that are activated during recovery from tobacco smoke exposure. Using whole genome gene-expression profiling of individuals initiating a smoking-cessation attempt, we sought to characterize the kinetics of gene-expression alterations in response to short-term smoking-cessation in the nasal epithelium. RNA was extracted from the nasal epithelial of active smokers at baseline and at 4, 8, 16, and 24-weeks after smoking-cessation and put onto Gene ST arrays. Gene-expression levels of 119 genes were associated with smoking-cessation (FDR < 0.05, FC ≥1.7) with a majority of the changes occurring by 8-weeks and a subset changing by 4-weeks. Genes down-regulated by 4- and 8-weeks post-smoking-cessation were involved in xenobiotic metabolism and anti-apoptotic functions respectively. These genes were enriched among genes previously found to be induced in smokers and following short-term in vitro exposure of airway epithelial cells to cigarette smoke (FDR < 0.05). Our findings suggest that the nasal epithelium can serve as a minimally-invasive tool to measure the reversible impact of smoking and broadly, may serve to assess the physiological impact of changes in smoking behavior.

Molecular subtyping reveals immune alterations associated with progression of bronchial premalignant lesions.

Bronchial premalignant lesions (PMLs) are precursors of lung squamous cell carcinoma, but have variable outcome, and we lack tools to identify and treat PMLs at risk for progression to cancer. Here we report the identification of four molecular subtypes of PMLs with distinct differences in epithelial and immune processes based on RNA-Seq profiling of endobronchial biopsies from high-risk smokers. The Proliferative subtype is enriched with bronchial dysplasia and exhibits up-regulation of metabolic and cell cycle pathways. A Proliferative subtype-associated gene signature identifies subjects with Proliferative PMLs from normal-appearing uninvolved large airway brushings with high specificity. In progressive/persistent Proliferative lesions expression of interferon signaling and antigen processing/presentation pathways decrease and immunofluorescence indicates a depletion of innate and adaptive immune cells compared with regressive lesions. Molecular biomarkers measured in PMLs or the uninvolved airway can enhance histopathological grading and suggest immunoprevention strategies for intercepting the progression of PMLs to lung cancer.